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ABSTRACT 


Ths  method  of  rtmdom  hazard  functions  Is  applied  to 
the  case  of  the  Welhull  distrihutioa  and  the  follow¬ 
ing  results  are  obtained:  It  Is  well  known  that,  In 
the  case  cf  a  Welbull  distribution  with  two  parameters 
a  and  8,  there  Is  no  way  of  estimating  and  testing  for 
the  shape  parameter  B  without  knowledge  of  the  scale 
parameter  a  using  the  usual  methods  based  on  maximum 
likelihood.  The  method  of  random  hazard  functions  Is 
now  used  to  obtain  a  consistent  and  asymptotically 
normal  class  of  estimates  for  8  Independent  of  any 
specification  whatsoever  on  a.  This  result  enables 
one  to  tes.  for  the  randomness  of  the  underlying 
failure  phenomena  under  the  Welbull  setup*  An  illus¬ 
trative  exaaq;ile  Is  given. 
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1.  IHTRODUCTION* 

The  main  object  in  reliability  analysis  is  to  estimate  for  each  specified  time, 
t  >  0,  the  probability  that  the  given  item  (where  item  denotes  a  component,  or 
a  subsystem  of  components,  or  a  system  of  subsystems)  survives  time  instant  t. 
This  function  of  t  is  usually  denoted  by  R(t).  The  con^lement  of  R(t),  which 
denotes  the  probability  that  the  item  failed  by  time  t,  is  denoted  by  F(t)  ■ 
l-R(t).  If  the  random  variable  T  denotes  the  observed  time  to  failure  of  an 
item,  F(t)  is  its  distribution  function,  which  is  also  referred  to  in  the 
literature  aa  "the  un'’*’*'>yine'  law  of  failures."  Assuming  that  F(t)  is  absolute¬ 
ly  continuous. 


dF(t)  -  f(t)dt  (1) 

there  f(t)dt  is  the  unconditional  probability  that  the  item  fails  during  the 
interval  (t,  t  *  dt).  The  function  f(t)  ia  called  the  probability  density  of 
the  underlying  lav  of  failures. 

Z(c)dt  denotes  the  conditional  probability  that  the  item  having  survived  time 
t  fails  during  the  next  interval  (t,  t  e  dt).  If  one  uses  this  notation,  the 
definition  of  conditional  probability  yields 

Z(t)dt-Yiy^  .  (?) 

Solving  Eqi,  (2)  tor  Tit)  gives 


t 

-J  Z(t)dt 
0 


r(t)  -  1-e 


(3) 
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Tte  function  Z(t}  It  called  the  hazard  rate.  The  (quantities  H(t),  Z(t),  and 
the  expected  value  of  the  random  variable  T  (representing  the  mean  time  to 
failure  [m.t.t.f.]  of  the  Item)  are  called  the  life  quality  or  item  effective- 
nett  parasietert.  It  should  be  said  at  this  stage  that  the  ultimate  end  In 
reliability  analytlt  it  to  estimate  or  establish  the  life  quality  or  Item 
effectiveness  parameters. 

Where  the  experimenter  has  sufficient  experience  vltL  the  Item  concerned  to 
specify  the  siathesiatlcal  form  of  the  underlying  lav  of  failures,  except  for 
the  knowledge  of  co^ain  parasieters  Involved,  the  statistical  estimation  pro¬ 
cedures  are  called  parametric.  Where  the  form  of  the  underlying  lav  of  failures 
cannot  be  aeevDaed  with  reasonable  accuracy,  thv  applicable  atatiatlcal  procedures 
tre  called  nonparametric. 

However,  there  are  many  significant  gaps  in  both  the  parametric  and  the  non- 
parsnitric  areas  of  reliability  estimation.  In  the  parametric  area,  exact 
methods  of  statistical  inference  are  pcssible  only  under  tne  exponential  lav 
of  failures.  With  respect  to  other  plausible  speclflcaticna  of  the  lav  of 
failures,  such  as  the  Welbull  diatributloo ,  method?  based  on  naximivn  likelihood 
end  order  etatlstice  provide  asymptotic  (large-s rumple)  statistical  procedures 
for  any  given  sampling  plan  (for  example,  random,  truncated,  and  censored). 

These  itethcds  are  usually  complicated  from  tim  point  of  vlev  of  evaluating  the 
•stloatcs  and  nuisance  parameters.  The  nonparanetrXc  approach  in  reliability 
is  relatively  nev  and  largely  unexplored. 

If  in  5q.  (2)  one  regards  Z(t)  not  es  a  function  of  the  real  variable  t 
but  as  a  function  of  tbe  tine  to  failure  T,  tt^en  one  lias  vhat  is  culled  a 
random  hassl'd  functiou.  At  this  nonect  ve  do  not  gc  into  detailed  properties 
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of  the  random  hazard  function  as  such.  However,  we  will  use  the  rendoo  hazard 
function  to  obtain  a  consistent  and  asymptotically  normal  class  of  estimators 
for  the  Weibull  shape  parameter  S  and  we  discover  that  the  sampling  properties 
of  our  estimate  for  B  are  asymptotically  free  from  a.  This  is  a  signlficanu 
property  which  the  class  of  estimators  for  B  given  by  the  method  of  maximum 
likelihood  does  not  (tnssess  even  for  large  samples. 

2.  ESTIMATION  OF  THE  WEIBULL  SHAPE  PARAMETER  B  BASED  ON  RANDOM  HAZARD  FUNCTION 
Consider  the  Weibull  distribution 

g 

F(tj  »  l-e““^  ,  a  >  0,  B  >  0.  (U) 

The  cumulative  hazard  rate  is  given  by 

t 

Y(t)  ■  j  2(x)dx  ■  at®  ,  (5) 

0 

Differentiating  (5)  with  respect  to  t  we  have  for  the  hazai*d  rate  Z(t)  of  the 
Weibull  distribution 


Z{t)  •  oBt®”^  ,  (6) 

Now  let  N  items  be  put  to  a  life  tasting  experiment  and  let  Tg,,..,  Tj^ 
denote  the  observed  times  to  failure  of  these  N  items.  Let  us  assume  that  the 
underlying  law  of  failures  is  Weibull  given  ty  (U).  Iq  other  words 
{1^}|  i*l,  2|.*.,N  are  independently  identically  distributed  random  variables 
with  conmoo  distribution  given  by  (b). 
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Bow  if  T  is  distributi^d  according  to  {k) 


'  o6  j 

0 

-  I  x^^'^e-^dx 
0 

«  a"^'^®r(l+K/6)  .  (7) 

where  F  (x)  denotes  the  gamma  function.  Consider  now  the  following  random 
hazard  function: 


g(T)  •  TZ(T)  ,  (8) 

where  T  is  the  observed  time  to  failure  and  Z(T)  is  the  random  hazaurd  function 
of  the  Weibull  dist^ribution  given  by  (6). 

Let  us  now  insider  the  following  random  variable  and  investigate  its  properties. 
Define 

N 

®'fZ 

i-1 

Before  we  go  any  further  let  us  note  that  g  is  not  a  statistic  as  we  ceuinot 

compute  its  value  as  soon  as  the  observations  are  available. 

We  have  „ 

N 

2  E(g(T^))  •  E(s(T))  £  (10) 

i-1 


U 


vbere  T  is  distributed  according  to  (U).  Nov 


E(g(T))  =  E(TZ(T)) 

»  EdaeT®"^)  (11) 

«  aeE(T®)  . 

a 

Substituting  for  ECT*^)  from  (7)  we  obtain 

E(g(T))  «  as  a“^  -  e  .  (12) 

A 

This  the  random  variable  S  ia  unbiased  for  S. 

A 

We  will  now  coiq;)ute  the  variance  of  the  random  variable  8.  From  (9)  ve  have 

Var(8)  -  I  Var(g  (T))  .  (13) 

Now 

E(g‘^(T))  -  E(aVT^®),  (ll») 

E(a^8^^®)  -  a^sV^rO) 

-  26^  .  (15) 

Combining  (12'^,  (13)  and  (15)  we  obtain 

Var(8)  ■  6^/N  (l6) 

A 

At  this  stage  we  have  discovered  that  the  random  variable  8  la  exactly'  unbiased 

2 

for  8  and  its  variance  8  does  not  dep«id  on  the  scale  parameter  a  of  the 

« 

VeibuU  distribution.  We  will  now  construct  a  statistic  8  with  the  following 
property,  viz; 

Pllra(6-S*)  -  0  ,  (17) 

N-*’  « 
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In  other  words  the  statistic  S  converges  in  probability  to  the  random  VEurlable 
0.  If  we  prove  (17)  t  proved  that  the  statistic  8  has  the  same 

asymptotic  properties |  in  fact«  has  the  same  as;/mptotic  distribution  as  the 
random  variable  6.  In  fact,  in  view  of  the  Lindeberg  and  Levy  version  of  the 
central  limit  theorem  <8ee  Cramdr  [l],  p.  213) 

A 

- ^H(0,1)  ,  (18) 

A 

in  other  words  6  is  asymptotically  normally  distributed  with  mean  0  and  variance 

2 .  » 

0  fS  independent  of  a.  Thus  if  we  can  define  a  statistic  8  suitably  and  show 

,  .  «  » 

(17)  for  that  6  then  we  havw  produced  a  statistic  8  with  all  the  desired 
properties  (asymptotically  or  for  Icurge  samples)  for  estimating  the  shape 
parameter  8  of  the  Welbull  distribution. 


Now  define 


where 


»'  ■  ?  S  ''i  • 

i»l 


V’'i>  ■ 


Rjj(Ti)  "I-  [number  of  olmervations  among  Tg,, . .  ,  (21) 

•  (22) 

and  is  a  sequence  of  nonnegative  constants  depending  on  the  sample  size  N 
such  that  Bjj  -►  »  as  N-*  •  and  -jp  0  as  N  and  finally  K(x)  is  a  window 

function  (see  Murthy  [2],  [3]  )  satisfying 


K(xi  S  0 
K(x)  »  KUx) 
lim  xK(x)  =  0 

Ixh  “ 

00 

J  i:(x)dx  =  1  , 


(23) 


Clearly  as  soon  as  one  chooses  the  sequence  {B^}  and  e  window  function  K(x) 
satisfying  the  above  properties  one  can  compute  the  statistic  g  given  by  (19) 
baaed  on  observed  values  of  time  to  failure  of  N  identical  items  put 

to  a  life  testing  experiment. 


Now  N 

■  I  S  ’'i  [''<■'!>  -V'l’l  • 

i»l 


It  is  evident  that  if  we  can  show 


(24) 


Plim  (Z(T  )  -  Z^(T  ))  -  0  ,  (25) 

00  ^  ^ 


it  follows  from  (24)  that 


Now 


Plim  (6-0*)  -  0 
N-^  • 


Vh> 

W  ’ 


(26) 


(27) 
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where  and  are  respectively  given  by  (22)  and  (21)  and 

^<v-^  . 

whore  f(t)  and  R(t)  are  respectively  the  density  and  reliability  functions. 

We  will  establish  that  converges  in  probability  to  Z(T^)  as  the  sample 

size  N  •  ly  establishing  that  the  denominator  of  converges  in 

probability  to  R(Tj^)  and  the  numerator  fg(Tj)  converges  in  probability  to  f(Tj^), 

We  will  prove  that 

Plim  (f  (T  )  -F(T.)  1  »  0  .  (29) 

N-*>  00  I  ”  ^  ^  ' 

Now  N 

■  I  •  <3°) 

where  U(x)  is  the  Heaviside  unit  function 

U(x)  ■  1  X  2:  0 

aO  other-rise  .  (31) 

We  have  i 

Fj^(Ti)  -F(T^)  -  I  j  1  ♦  U(T^-T^)  ♦...+  U(Tj-T^_^)  ♦ 

+...♦  U(T^-Tjj)  -NF(Tj)  f  .  (32) 

Since  F(T^)  is  uniformly  distributed  in  ((^,1)  ' 

E(F(T^))-^  .  (33) 
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Also  since  and,  say,  are  independently  identicedly  distributed 

E(U(T,.Tj)  «  P(T,  s  Tj  (3U) 

11  i  X 

1 

'  2  • 

Using  (33)  and  (3^)  and  taking  expectation  on  both  sides  of  (32)  we  discover 
that 

1+  i(N-l-N) 

e(  I-„(Tj)-F(Tj)|  .  - 5 - 

1 

“  2N  (35) 

and  hence 

lim  E  I  Fjj(T^)-F(T^)j  »  0  .  (36) 

Now 

Fjj(T^)  -  F(T^)  -  i  [u(T^-T^)-F(T^)]  +...+  [u{T^-T^_^)-F(T^)  ] 

+  [u(T^-T^^^)-F(Tj)j  ♦...+  [u(T^-Tjj)-F(Tj)] 

+  1  -  F(Tj) 

(37) 
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Write 

-  U(T^-Tj^)-F(T^)  , 

Zg  -  U(T^~T2)-F(T^)  , 

•  •  • 

^1-1  “  , 

h  -  U(Ti-Ti^i)-F(T^)  , 


(36) 


h.i  “  • 


Then 


Bw 


Var  [Fj,(Ti)-F(T^)  ]  «  jA  [variZj^+Zg  ♦...+  Zj^^^  ♦  l-F(T^) }]  {39) 

“7  Var(Z^  ♦...♦  Zj^_l) 

*Cov^(,i.r(Tj),Zj*Zj*....  z^_\ 

Var(Z^)  -  V«r(U(T^-T^))  ♦  ^  -2Cov  |  F(Tj) .UfT^-T^)]  (Uo) 

.  1  ♦-!  1  .  1 

u  12  6  ?r  • 

•  •  • 

Car(Z^.Zj)  .  J  I  (« '  |  W 

:<  <lF(T^)dr{T2)dF(T^)  , 


COO 
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since  E(Zj^)  »  0,  i  »  1, 
Expanding  (Ul)  we  obtain 


CovCZ^.Zg)  “  [  j  j 
0  0  0 


U(  )  U(  T^-Tg )  dF(  )  dF(  Tg )  dF(  ) 


-IIJ  F(  )  U(  T^-Tg )  dF(  )  dF(  Tg )  dF(  ) 

0  0  0 

F(T^)U(T^-Tj^)dF(Tj^)dF(T2)dF(T^) 

0  0  0 

00  a» 

•H!  F^(T^)dF{Tj^)dF(T2)dF{T, ) 

0  0  0 

-  l-E{F(Max(T^,T2)))  -  E(F(nax(T^,Tg) )) 

Now 

m  m 

E(F(max(T^,T2)l)  -  J  r(x)2F(x)f(x)dx  -  2  J  F^{x)l'{x)dx  - 

0  0 


(1*2) 


(1*3) 


i: 


Combining  {k^.)  and  (I43)  ve  discover  that 


Cov(Zj^,Z2)  ■  0 

Similarly 

Var(Zj^)  '  ^  ,  i  “  1,2,...  ,11-1, 

Cov(Z^,Zj)  -0  I  j<  J 

Finally 

Cov(Z^,(l-F(T^))  -  E(Z^(1-F(T^))  [since  E(Z^)  -  0  ] 

-  E(  [l-F(T^)]  [u(Tj^-Tj^)-F(T^)] ) 

-  E(U(T^-T^)-F(T^)U{T^-T^)-nT^)+F^(T^)) 

-i.i.i.i.O  . 

2323 


m 


(U5) 


(j>^) 


Siffillcurly 

CoY[Zj,(l-F(Tj))j  •  0  ,  i  m  1,2,...,N-1  . 

Cooblning  (39),  (*»5)  and  {kf)  we  discover  that 


(U7) 


(J*8) 
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Hence 


lim  Var  [fjj(T^)-F(T^)] 


0  . 


(»49) 


Conbinlng  (36)  and  (U9)  ve  finally  discover 

Plim  (F„(T,)-F(T,))  -  0  . 

N-*'  « 


which  Isqplies 


Plim  (R_(T.)-R(T,))  -  0 


(50) 


To  conqplete  the  proof  of  (25)  we  now  have  only  to  show  that  fj|{Tj^J  converges 
in  probability  to  f(Tj^), 


In  other  words  we  have  to  show  that 


Plim  (f_(T.)-f(T.))  ■  0 


(51) 


where 


Let 


J-1 

•  K(yT^-Tj))  . 


(52) 


Vl  -  *  <53) 

9  9  9 
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Then  f„(T. )  cea*  be  written  as 
N  1 


B  K(0)  B 

f  (T  )  ■  -JL— —  ♦  >  ^4 

J-T- 


Coosider  norj 


E(fjj(Tj)-f(T^)) 


B  ■ 


N-1 


B^K(0)-1(T  )  B 

J-1 


B  E(f(r))  3„ 


N-1 


Bincfa  by  aasumptioa 


i^o) - ji-.  IS  «V!T,)) 

B.,  E{f{T  )) 

-UK(0) - 


®R 

— -  0  M  N  -  •  and  K{0)  is  a  constant, 


lim 

R  - 


K{0)Bjj 


0  , 


Now 


|E(f(Tj))  -  I  J  f(t)f{t)dt 


f^(t)dt  . 

0 


Ik 


(5k) 


(55) 


(56) 


I 


Assiuning  that  » 

[  f^(t)dt  <  -  , 

0 

which  is  certainly  true  in  the  ease  of  the  Weibull  diatribution,  we 


lim  i  E{f(Tj))  -  0  . 


Now 


E(Bjj(Xj^«.f(T^)))  ■  J  j  Bj^[K(Bjj(T^-T^))-f(T^)ldF{Ti)dF(Tj) 


0  0 


j*  dF(T,)  I  a^[K(Bjj(T.-T^))-f(T^)]f(Tj^)dTj 


It  now  follows  from  a  leona  (see  Murthy  [3],  p.  1028)  that 

m 

lim  j  B^K(Bj,(T^-Tj))-f(Tjn  f(T^)dT2^  -  0  . 

0 


Coahining  (59)  ard  (60)  we  obtain 

11«  E(Bj^(X^-f(Tj))  •  0  . 


(57) 

obtain 

(58) 


(59) 


(60) 


(61) 
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Finally  combining  (55) »  (56),  (58)  and  (6l)  we  discover  that 

lira  E(fjj(T.)-f(Tj^))  =  0  .  (62) 

N 

Appealing  to  the  lemma  (see  Murthy  [3],  p.  1028)  repeatedly  and  carrying  out 
a  straightforward  but  laborious  computation,  one  discovers  that 

lim  Var(fj^(Tj^)-f(T^))  *  0  .  (63) 

N  • 

Combining  (62)  and  (63)  we  obtain 

Plim  (fjj(Tj^)-f{T^))  «  0  .  (6U) 

N  -*•  * 

In  view  of  (50),  (64)  and  (25)  we  finally  prove  that 

Plim  (6-6*)  «  0  .  (65) 

N  • 


Hence 

THEOREM  « 

The  sequence  of  statistics  converges  in  distribution  to  a  normal 

P 

distribution  with  zero  mean  and  unit  variance  as  the  sample  size  N  tends  to 
infinity. 
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3,  TFST  FOR  EXPONEBTIALITY  IH  THE  WEIBULL  SETUP 

The  Weibull  distribution  given  by  (U)  reduces  to  the  exponential 
distribution  if  the  shape  parameter  8  1. 

To  teat 

Hq;  6=1  ,  (66) 

we  use  the  theorem  Just  established  as  follows.  Suppose  that  the 
sample  size  H  is  large  ei.ough  for  the  normal  approximstion  proved  in 
the  theorem  to  be  valid.  Choose  a  sequence  {Bj.}  and  a  window  K(x). 
Compute  6  and  (6  -l)  If  the  latter  value  exceeds  the  QJ?  value  of 

a  normal  distribution  with  zero  mean  and  unit  variance  reject  the 
hypothesea  Hq. 


U.  M  EXAMPLE 

Consider  the  estimate  6*  of  6  given  by  (19) »  namely; 


“  1 1 

i=l 


(67) 


where  T,  ,  T_,...,T„  are  the  observed  times  to  failure  and 

1’  2’  *  R 


(68) 


njj(Ti)  =  ^  (number  of  observations  among  T^,  >  Tj^j  , 

(69) 

N 


J-1 


(70) 


17 


emd  is  a  sec^uence  of  noonegative  constants  depending  on  the  sample  size 
N  such  that 


lim  “TT  »  0 

K  ott 


(71) 


and  finally  the  window  K(x}  satisfies  (23). 


To  illustrate  the  use  of  the  statistic  to  test  the  hypothesis  given  by 
(66),  let  us  consider 

Bf,  =  log  N  (72) 

K  (x)  »  i  |x|s  1  , 

^  * 

«  0  otherwise; 


clearly  given  by  (72)  satisfied  (71)  and  K(x)  defined  by  (73)  satisfies 
the  window  condition  (23). 


If  we  substitute  (72)  and  (73)  in  (67)  the  statistic  g  in  this  illustration 
reduces  to 


1+  I  number  of  observations  Tj(j  ^  l),  such  that 


*  _  lOfi 


11*1 


i-1 


number  of  observations  T 


J 


such  that  T 


j '  *1 1 
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As  soon  as  the  observations  are  avedlable  the  statistic  (7^+) 

can  at  once  be  computed  and  the  procedure  for  rejection  or  acceptance  of  the 
hypothesis  can  be  carried  out. 
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